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We present a comprehensive, nonperturbative analytical method to investigate the dynamics 
of time-dependent oscillating scalar field configurations. The method is applied to oscillons in a 
0^ Klein-Gordon model in two and three spatial dimensions, yielding high accuracy results in the 
characterization of all aspects of the complex oscillon dynamics. In particular, we show how oscillons 
can be interpreted as long-lived perturbations about an attractor in field configuration space. By 
investigating their radiation rate as they approach the attractor, we obtain an accurate estimate of 
their lifetimes in d = 3 and explain why they seem to be perturbatively stable in d = 2, where d is 
the number of spatial dimensions. 



I. INTRODUCTION 

Nonlinear field theories contain a large number of local- 
ized solutions that display a rich array of properties [T] . 
Of particular interest are those that are static and sta- 
ble, that is, that retain their spatial profile as they move 
across space or scatter with each other, as is the case of 
sine-Gordon solitons. The details of the solitonic config- 
urations are, of course, sensitive to the dimensionality of 
space and to the nature of the field interactions. Long 
ago. Derrick has shown that, in the case of models with 
just a real scalar field, no static solitonic configurations 
can exist in more than one spatial dimension [2]. Given 
that most models of interest in high energy physics in- 
volve more complicated fields in three spatial dimensions, 
this restriction was somewhat frustrating. Fortunately, 
the subsequent exploration of a variety of models led to a 
plethora of static, nonperturbative, localized field config- 
urations. Examples include topological defects, solutions 
of models usually involving gauge fields that owe their 
stability to the nontrivial topology of the vacuum, such 
as strings and monopoles 'W , and the so-called nontopo- 
logical solitons, solutions of models where a conserved 
global charge is trapped inside a finite region of space 
due to a mass gap condition, such as Q-balls |4] and the 
models with a real and a complex scalar field of Fried- 
berg, Lee, and Sirlin 0. 

In the mid-nineties [S] , a new class of localized nonper- 
turbative solution began to be explored in detail, after 
being proposed earlier [7]. Named oscillons, such long- 
lived solutions have the distinctive and counter-intuitive 
feature of being time-dependent. In spite of this, the 
nonlinear interactions act to preserve the localization of 
the energy, which remains approximately constant for a 
surprisingly long time [8J. During the past few years, 
oscillons have attracted much interest. Their properties 
were explored in two [9] and higher [10^ spatial dimen- 
sions, in the presence of gauge fields [TT], in the standard 
model of particle physics [T2], and in a simple cosmolog- 



ical setting [TF . There have also been detailed attempts 
at understanding some properties of oscillon-related con- 
figurations (typically with small-amplitude oscillations), 
including their longevity, using perturbative techniques 
|14j . On the other hand, a treatment explaining the re- 
markable longevity of oscillons in models related to spon- 
taneous symmetry breaking, and thus of obvious interest 
in particle physics and cosmology, has been lacking. The 
situation was partially remedied recently, when we pub- 
lished a preliminary treatment of the problem |15j . In 
the present work, we greatly extend the range of our dy- 
namical theory of oscillons in scalar field models, include 
the details of many key derivations and demonstrate its 
accuracy in reproducing numerical results. Our approach 
is general enough to be extended to different scalar field 
models that exhibit long-lived, time-dependent localized 
configurations. 



II. LINEAR VS. NONLINEAR DYNAMICS AND 
THE OSCILLON MASS GAP 

In order to introduce some of the basic quantities 
needed for our theory, it is instructive to start by review- 
ing some of the main properties of relativistic oscillons. 
We will do so in the context of a simple (j)* model with a 
symmetric double-well potential, as this is also the main 
focus of the present work. To begin, consider the La- 
grangian for a spherically-symmetric, real scalar field in 
d-spatial dimensions. 
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where V{(j)) = rn^cj)^, and q = 27r'^/2/r(d/2) is the unit- 
sphere volume in d dimensions. Quantities are scaled 
to be dimensionless as follows: (j) = m^''~"'^^/^0o and 
r'' — rQ/m. We will henceforth only use dimensionless 
variables, dropping the subscript "0" . 

We start by investigating the linear theory so that we 
can more easily contrast it with nonlinear models that 
give rise to oscillons. Since oscillons have been shown 
to maintain their approximate Gaussian-shaped spatial 
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profiles during their lifetimes, we will write the scalar 
field as 



{r,t) = A{t)P{r;R) = A{t)e 



(2) 



Here, is the time-dependent amplitude of the con- 
figuration and P{r\ R) its spatial profile, which is param- 
eterized by the radial extension R. 



A. Linear Dynamics 

As shown in Rcf. [8J, the d = 3 linear theory with a 
Gaussian-profile initial condition has the solution. 
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where A is an arbitrary initial amplitude and the disper- 
sion relation is uj ~ (fc^ + ^jf^-^^ss)^ i where the mass fre- 
quency Wmass = %/2. To Calculate Fun, the decay width 
associated with the above solution, we recall that in ^8. 
it was also shown that Eq. [3] can be approximately inte- 
grated to obtain (at r = 0, the configuration's maximum 
amplitude) 
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whose envelope of oscillation is given by (/)(0, t) = Aq/{1+ 
which reaches 1/e of its initial value in a time 
given by Tiincar — ■^'i&i^ma.ssR^ ■ This yields the linear 
decay width Fiin, 
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In the linear theory, any initial configuration-or exci- 
tation above the vacuum-will quickly decay by emitting 
radiation. The key difference between the linear and the 
oscillon-supporting nonlinear models is that, in the lat- 
ter case, the decay modes are strongly suppressed. It 
is this suppression that gives rise to the oscillon's re- 
markable longevity. Our goal in this paper is to make 
this statement quantitatively precise. To obtain the lin- 
ear radiation distribution, which we denote by we 
simply take the fc-space representation ke^^ 1^ of the 
Gaussian in Eq. |3]and express it in terms of uj using the 
dispersion relation uj — (fc^ -1-2)5: 



6(cij) = ^[cjje 



(6) 



hiiS) is a lopsided distribution with frequencies above 
^mass = \/2 and peaked at w^ax = (2 + 2/i?^)^/^. Now 
define uJidt and Wj-jght to be the frequencies where the dis- 
tribution h(uS) rises to half of its peak value. A straight- 
forward calculation gives wicft ~ (2 -I- .203/i?^)2 and 
Wright — (2 -I- 7.38/i?^)2. The radiation distribution is 
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FIG. 1: Schematic description of the hnear radiation peak 
centered on ujin being shifted to the left by the presence of 
nonlinearities. An oscillon will form if the peak is shifted 
far enough into the nonlinear region such that it no longer 
overlaps significantly with the linear peak. 



then approximately centered on the frequency wnn, given 
by (see Fig. [l|. 



(^Icft + bright) 

.203\^ 1 



(7) 



7.38 



which we take to represent the dominant linear radiation 
frequency. 



B. Nonlinear Dynamics and Decay Rate 

Imagine now that one or more nonlinear terms are 
added to the linear potential and that the field is again 
initialized with the same localized Gaussian perturba- 
tion. Roughly speaking, the nonlinearities will shift the 
dominant linear oscillation frequency tj^nn, to a new value, 
denoted by io^\. If the added terms serve to decrease the 
curvature of the potential, then Wni < (^lin- Such a situa- 
tion is depicted qualitatively in Fig. [T] where the arrow 
indicates how the shift in frequencies occurs. 

When nonlinearities are efficient enough that is 
lowered substantially below o^un, an oscillon may form. 
The reason for this is the following. As the initial field 
configuration begins to oscillate, it will attempt to emit 
small-amplitude radiation waves in an effort to dissipate 
its energy. However, if is sufficiently less than wun, the 
bulk of the frequency components composing the oscil- 
lation will be unable to excite small-amplitude radiation 
waves (since the configuration can only radiate appre- 
ciably in the frequency range cjictt < < Wright)- This 
condition leads to the stabilization mechanism responsi- 
ble for the formation of oscillons. 

It then follows that the oscillon will enter the nonlin- 
ear regime if the two peaks in Fig. [l]do not significantly 
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overlap. Mathematically, the right "edge" of the nonlin- 



ear peak, given by Wni 
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must be less than the left 



"edge" of the linear peak, given by cjun — ^Tu^- There- 
fore, we have uini + iTni < i^iin — -^^Un which, by defining 
(jjgap = wiin — Wni as the frequency gap between the linear 
and nonlinear peaks, becomes 
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> - (r„i + Tiin) 
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Since nonlinearities tend to increase a configuration's life- 
time and thus decrease its decay rate, it follows that 
< Fni < Fiin. One can thus state that if 



^gap 



> r 



lin ; 



(9) 



then the configuration will be forced into the nonlinear 
regime. In other words, Eq. [9] is a necessary condition 
for the formation of an oscillon. As we will soon see, if 
an oscillon is formed and, during the course of its time 
evolution, reaches a point where Eq. [9] is no longer sat- 
isfied, it will cease to exist. This implies that oscillons 
decay when 
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To obtain the nonlinear radiation frequency uj^i and 
thus ti-'gap, we substitute Eq. [2] into Eq. ^and integrate, 
giving 



2 ^(^max) 



(11) 



where now includes nonlinear terms and E is the 

energy which is found by taking the appropriate Legen- 
dre transform of the Lagrangian and evaluating it at the 
upper turning point of an oscillation, ^max- The oscil- 
lation frequency of the oscillon, Wni, at a given time is 
given by 
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dt = 2 
Ja^^^ a 



(12) 



where A = [2E/cr - 2V{A)]^/^, cr = {n/2Y/^R'^, and 
is given by V{Ajnin) = ^(^max)- 



C. The Attractor Point 

Results from numerical simulations suggest that there 
exists an attractor point in configuration space to which 
the oscillon tends. It was noted in [TS] that one can ob- 
tain the energy of this attractor point (in cf)'^ models) by 
finding the minimum energy which has the property that 
the effective potential V{A) possesses at least one point 
for which V"{A) < 0. In order to compute the attractor 
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FIG. 2: Minimum oscillon energy as a function of core am- 
plitude for a double- well potential in d = 3. The minimum of 
this curve is the attractor point with Eoo — 37.69. 



point, it is easier to work within a specific model. Choos- 
ing V{A) — (j)'^ — (j)^ + (/)'^/4, we obtain, using Eq. [2] and 
integrating over all space. 
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F"(^).(2+A)-6(?)^ + 3|^ (14) 

Equate this to zero and solve for i? as a function of A. 
Then substitute the result into Eq. [TT] to eliminate R, 
yielding energy as a function only of A. This curve pos- 
sesses a minimum, shown in Fig. [2] for d = 3, the energy 
of which yields the correct attractor energy Eoo of the 
oscillon, which has numerical values E^o — 4.44 in d = 2 
and Eoo — 37.69 in d = 3. 

Given the energies calculated above, there is a locus of 
points in {A, R) parameter space (the thinner solid lines 
in Figs. [3] and [4]) which possess these asymptotic values 
for the energy; one of them is the attractor. To locate 
the attractor point, we need to determine its amplitude 
coordinate, denoted A^o- In d = 2, this is most easily 
determined numerically to be — -3. In d = 3, we 
cannot determine A^o numerically since the oscillon de- 
cays before reaching it. Therefore, in d = 3, we must 
estimate A^o analytically. To do this, we choose the 
point (satisfying E — Eoo) which has Wni — t^^mass, even 
though, in reality, its frequency is slightly less than Wmass 
(never above it). As we will see, this approximation will 
suffice for our purposes. From Eq. |12[ the amplitude 
which gives ojni — '^-'mass (in d = 3) has numerical value 
Aoo ^ .456. 

Given the pair {Aoo, Eoo) — (.456, 37.69) in d = 3 and 
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(.3,4.44) in d = 
. We obtain Ro 



we can use 
4.79 for d = 



Eq. 

= 3, 



and 



4 




i?oo ^ 5.77 in d = 2. The circles in Figs, [s] (d = 3) and 
|4] (d = 2) mark the locations of the attractor points. 
Given Eqs. Is] [tI and 12 (and their equivalents in d = 2 



Amplitude 

FIG. 3: The thick solid hne represents the locus of points 
satisfying the condition LOga,p = Tun in d = 2 ("line of exis- 
tence") calculated analytically. The thin solid line represents 
the locus of points that have the attractor energy E ~ 4.44. 
The dashed line is the numerical minimum radius based on 
Gaussian initial configurations. Oscillons can only exist if 
their energy E > Eao and if they have a core amphtude and 
average radius lying above the line of existence. This is why 
the numerically measured minimum radius follows the line of 
existence for A > 1 but follows the line of minimum energy 
for A < 1. The arrow indicates the value of the (constant) 
minimum radius calculated in [8]. The circle represents the 
location of the attractor point; since it lies above the line of 
existence, oscillons will be absolutely stable in this system. 
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FIG. 4: The thick solid line represents the locus of points 
satisfying the condition ojgap = Fiin in d = 3 ("line of exis- 
tence") calculated analytically. The thin solid line represents 
the locus of points that have the attractor energy E ~ 37.69. 
The dashed line is the numerical minimum radius based on 
Gaussian initial configurations. The arrow indicates the value 
of the (constant) minimum radius calculated in 8_. Oscillons 
can only exist if their core amplitude and average radius are 
above the curve (wherein E > Eao is automatically satisfied). 
The circle represents the location of the attractor point; since 
it lies below the line of existence, oscillons will not be abso- 
lutely stable in this system. 



shown in Appendix A) we can also calculate the quan- 
tities in Eq. |9] that is, the oscillon existence condition, 
as a function of the parameter pair {A,R). In Figs, [s] 
(d = 2) and|4] [d — 3), the thicker continuous lines rep- 
resent the locus of points satisfying Wgap — defining 
the boundary line between the region where oscillons may 
exist (above the line, and provided that E > Eao) and 
where they cannot exist (below the line). We will often 
refer to this boundary as the "line of existence." As a test 
of the existence condition (Eq. [9]), we also plot the nu- 
merical result for the "minimum radius" as a function of 
amplitude (dashed line), found by pinpointing the min- 
imum initial radius which causes a configuration to live 
longer than the linear decay time. The arrows in Figs. [3] 
and|4]indicate the values of the minimum radii calculated 
in [8J which possess no amplitude dependence and thus 
provide only limited information. 

As can be seen, the attractor point in d = 3 lies below 
the line of existence curve, explaining the finite lifetimes 
of oscillons in that system: the configurations decay be- 
fore reaching the attractor point. On the other hand, 
the attractor point in d = 2 lies above the curve, explain- 
ing the seemingly infinite oscillon lifetimes observed in 
numerical simulations. We can thus interpret the oscil- 
lons as time-dependent perturbations about the attractor 
point. Those in d = 3 are unstable, albeit some can be 
extremely long-lived. Those in d = 2 are at least pertur- 
batively stable. 

In situations such as d = 3, where oscillons eventu- 
ally decay, it is interesting to compute their lifetimes and 
how they depend on their radiation rate. In a recent 
work, we presented the basic features of a method de- 
signed to do so [T5]. In the next section, we develop 
the appropriate formalism in detail, based on the over- 
lap between the nonlinear and linear radiation spectra. 
We point out that our formalism is, in principle, appli- 
cable to any time-dependent scalar field configuration, 
offering a much-needed handle on how to compute radia- 
tion rates of nonperturbative configurations in relativistic 
scalar filed theories. 



III. LIFETIME OF LONG-LIVED OSCILLONS: 
GENERAL THEORY 

In the situation that ujgap > and an oscillon has 
formed (above the solid curve in Fig. |4]) it will begin 
to radiate small amounts of energy. In this section, we 
will derive a general equation governing its radiation rate 



so that, in the event that 



— Fun and the oscillon 



'^gap ^ lii 

decays, we may calculate its lifetime. As in the previ- 
ous section, our general approach will be to compute the 
overlap between the nonlinear peak and the linear peak. 
Since we are assuming that Wgap > ^ T^i for a long- 
lived oscillon, this overlap will be small; the amount by 
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Frequency 

FIG. 5: Schematic of an oscillon frequency distribution show- 
ing the tail penetrating the radiation region. The graph is 
plotted in units of Wmass. In the inset we show a close-up 
view of J-{uj) and respectively (dashed lines), and their 

product Q.{uj) (solid line) for the system given by Eq. [l] in 
d = 3 for typical values of the various parameters. Note that 
the curves have been vertically scaled so that all are visible 
on the same graph. 



which it differs from zero will determine the radiation 
rate. 

In Fig. |5]we plot a possible frequency distribution for 
an oscillon (the width of the peak centered on Wni has 
been greatly exaggerated for visibility) in units of the 
mass frequency (wmass — V^)- Note how the tail of the 
distribution "leaks" beyond the mass frequency. It is this 
leakage that will determine the radiation-rate and thus 
the decay rate of the oscillon. 



A. The Long-Lived Oscillon Radiation Equation 

In order to compute the oscillon decay rate, we model 
oscillons as spherically-symmetric objects whose radia- 
tion obeys a distribution (amplitude per unit frequency) 
n(a;) which consists of a narrow peak of width 5 centered 
at some frequency cji-ad- The radiation flux <i> (energy per 
unit time per unit surface area) emitted by such an object 
is 
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where S is the surface area of the oscillon, p is the 
radiation- wave energy density, v — Wrad/^rad is the phase 
velocity of the wave, fc^ad is the wave number, and A is 
the amplitude of the radiation wave, which is given by 



A = 



^(LS)(ko ~ ri(a;i.ad)<5. 



(16) 



The function ^{uj), which represents the amplitude per 
unit frequency of the radiation wave, is simply deter- 
mined by the "overlap" between the oscillon and the lin- 
ear radiation peaks. Taking ^(uS) (the overlap function) 
to be the product of the nonlinear peak and the linear 
radiation distribution b(uj) obtained in the previous sec- 
tion, we have 



T{uj)h{uj) 



(17) 



where .F(tj) is the nonlinear peak (Fourier transform of 
the oscillon 's core), a is a proportionality constant to be 
determined, and b = ab{uj). The inset in Fig. [s] shows 
a typical overlap function for an oscillon m d — i. (See 
Appendix C for details.) In Appendix B we show that, 
in the tail. 
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Combining Eqs. [15] [17] and [18] and letting S 
CdR'^^^, we have 



dE fdAy 
where is a time-dependent parameter given by 
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(19) 



(20) 



Eq. [19] is a differential equation which must be satisfied 
by a long-lived oscillon. In Appendix C we compute both 
a and S in general and in the context of the model of Eq. 

m 



B. Integration of the Long-Lived Oscillon 
Radiation Equation 

In this section, we will attempt to integrate Eq. [19] to 
obtain the oscillon energy as a function of time. This pro- 
cess is not straightforward since Eq. 19 contains deriva- 



tives of two different quantities (amplitude and energy) 
and the parameter r] possesses a complicated time depen- 
dence which is not known. However, in Appendix D we 
develop a simple method to solve this problem, based on 
the assumption that the timescale associated with the os- 
cillon's loss of energy is closely related to the timescales 
associated with the rates of change of all other oscillon 
parameters. The result of our approach is Eq. D4 which 
we will employ below. 

First, write Eq. |D4| for the cases X — A and X = rj 
and substitute into Eq. [19] obtaining. 
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Now, consider the situation where r/oo — (which is the 
case in the system we are studying here since, in d = 
3, the attractor point satisfies lu^^i ~ Wmass, leading to 
(jjrad — Wmass, which causes T] to bc zero there). In this 
situation, Eq. 21 is somewhat simphfied: 



where the constant 7^ is given by 

7r, 



7ij 
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and the exponent is (using that [E]oc = 0) 

= 2{1~ pa) + Pr,. (24) 



Eq. 22 is a constant-coefficient, ordinary differential 
equation governing the oscillon energy as a function of 
time, and can be easily integrated: 



Eit) = E, 



Ei — En 



[l + ^f,g{E,-Eno)n]' 



(25) 



where g = p^ — 1 and Ei is the energy at t = 0. Eq. 25 
is the energy of an oscillon as a function of time. 

As shown in Fig. |4] in d = 3 an oscillon will always 
decay before reaching E — E^o- The decay is quite sud- 
den, a burst of scalar radiation. As stated in section II, 



this occurs when w, 
denoted i?D, is given by 



gap 



Fiin; hence the decay energy. 
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To calculate the lifetime, denoted Twic, which is defined 
as the amount of time taken for the oscillon to decay from 
a sufficiently high initial energy down to Ed, we invert 
Eq. |25]to yield time as a function of energy and evaluate 
at Ej^: 
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When Ei — E^o ~> E — E^q, the function t{E) tends to a 
finite, maximum value leading to 
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(28) 



This means that, when the initial energy Ei is much 
larger than the energy £"13 in question, the time it takes 
for E{t) to fall from Ei to £'d becomes independent of the 
initial condition (i.e., from Ei). One can then say that, 
in a restricted sense, the long-lived oscillon is decoupled 
from initial conditions: if the necessary conditions for its 
existence are satisfied, a variety of initial configurations 
will approach an oscillon. Recent studies that have ob- 
served the emergence of oscillons from stochastic initial 
conditions after a fast quench offer strong support for 
this claim fTO. 
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1 (top curve) and 
plotted against the refer- 



FIG. 6: Analytical calculation of g = 
Pa (bottom curve) using Eq. 
ence amplitude Ar along the line of existence. This gives the 
spectrum of possible values of g and pA across the various os- 
cillons in this system. The dots mark the theoretical values of 
g ~ 2.67 and pA — .66 assumed by the longest-lived oscillon, 
which has Ar — .92 (thicker line in Fig. [8|. 



Eqs. [26| and [28] together give the lifetime of an oscillon 
and can be considered the main results of this paper. 
Before moving on, we note that, for a long-lived oscillon, 
F„i is given by 



Tni = r^i 
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where we've used Eqs. [Bg) |D3] and|22| This is related 
to the lifetime by (combine Eqs. 28 and 29 1 
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C. Sample Calculation: Klein-Gordon Field in 

d = 3 

We now apply the above results to the system given 
by the Lagrangian in Eq. [T] for d = 3, supplemented by 
the nonlinear potential of Eq. [T3j We will begin with 
the existence condition shown in Fig. |4] In section II. c 
we calculated the coordinates of the attractor point and 
obtained (^00,^00) — (.456,4.79). Comparison with Fig. 
|4] reveals that the attractor point lies below the curve, 
and thus that Wgap/Fjin < 1 there. As we noted before, 
stable oscillons will not exist in this model. However, 
there are still long-lived oscillons obtained by initializing 
the field sufficiently far from the attractor point. 

As these structures radiate energy, their amplitude A 
and radius R will change in time. Hence, they will trace 
out trajectories in the {A, R) plane of Fig. |4] all of which 
will eventually intersect the line of existence. To calcu- 
late these trajectories, begin by writing Eq. |D4|for A{t) 
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FIG. 7: Analytical ca lcul ation of 7^ (top graph) and 7a (bot- 
tom graph) using Eq. 33 plotted against the reference ampli- 
tude Ar along the line of existence. This gives the spectrum 
of possible values of 7^ and 7a across the various oscillons in 
this system. The dots serve to mark the theoretical values of 
7^ ~ 2.0 X 10~® and 7a — -21 assumed by the longest-lived 
oscillon, which has Ar — .92 (thicker line in Fig. [8|. 



and R{t), obtaining [A{t) - A^] = -/A[E{t) - E'oo]^'* and 
[R{t) - i?oo] = 7R[E{t) - E^]P'^, respectively. Substitut- 
ing the first into the second to ehminate the energy, we 
obtain 



R{t) = Roo+ lB.A[A{t) - Aoo]"^-' 



(31) 



IbJYa^ and pba 



PrIpa- Note that 



where 7_ra = ihi ia 
when A — Aoo , R — Roo ■ 

Since each possible trajectory will intersect the line 
of existence, every point along the line of existence is 
a point along some trajectory. Thus, we can choose an 
arbitrary point on this line (which we will call the "refer- 
ence point") and use its coordinates, labeled {Ar, Rr), to 
calculate all of the dynamical exponents {pAi PR, P^^ Pu>, 
Pr^, etc.) associated with the trajectory that intersects 
that point. Eq. |D12| reproduced below, was evaluated 
numerically for X ~ and X — pA, respectively: 



Px 



\ OVe ) OVe \ OVA 



dvx 

dvA 



(^ 



duE ) 



(32) 



Fig. |6] shows the results for the exponents g = p^ — 1 
(top curve) and pA (bottom curve) as a function of Ar- 
Given values for px , we can calculate by evaluating 
each side of Eq. D4 at the reference point (wherein X 



assumes the value X^ ~ X[Ar,Rr\) and solving for 7^, 
obtaining 



Xr — Xn 
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Fig. [7] shows the results of such a calculation for 7^ (top 
graph) and "fA (bottom graph). 



FIG. 8: Various oscillon trajectories. Note that they all 
tend to the attractor point but intersect the line of existence 
(dashed line) before doing so. The thicker trajectory marks 
the longest-lived oscillon (it intersects the line of existence at 
Ar ~ .92 and Rr ~ 3.25). 



Having computed px and 7^ for the parameters A 
and R allows us to calculate the coefficients in Eq. [3llfor 
each trajectory (i.e 



32 and 33 for various 
The result 



compute Eqs. 
reference points along the line of existence) 
of this is shown in Fig. [8] for selected trajectories. Note 
that they all asymptotically tend toward the attractor 
point but intersect the line of existence before doing so. 

As an example, consider the trajectory which intersects 
the line of existence at Ar ~ .92 (this will be shown to 
correspond to the longest-lived oscillon in the system of 
Eq. [1]). From Figs. [6] and [t] we have g ~ 2.67 and 7^ ~ 
2.0 X 10^^. If this oscillon were initiated at an energy 
of Ei ~ 82.5 (which will be the case for the numerical 
simulation we will be comparing to) then Eq. 25 for that 
oscillon becomes 



E{t) ~ 37.69 



44.8 



[1 + (.136)i] 



.375 ' 



which, for times t > 100 is approximately 

94 77 

E{t)-E, 



i-375 



(34) 



(35) 



Since the decay energy for this oscillon is E\) = Er = 
Er{Ar,Rr) ^ 41.0963, Eq. [28]yields 



1 



fo 



[2.0 X 10-6] [2.67] [41.10 - 37.69]2-67 
7100. 



(36) 



Given Eq. |34]we can write, for example, an expression 
for the oscillon amplitude as a function of time. From 
Figs. [6]and|7]we have pA — .67 and 7^ — -21. Combining 
this information with [A — ^00] = 1a[E — Eaa]'''^ and Eq. 
IMlwe have 



A{t) ~ .456 



2.58 
[1 + (.136)^] 



.249 ■ 



(37) 
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FIG. 9: Analytical results for critical values of the ampli- 
tude (top) and radius (bottom) along the line of existence 
(dashed lines) as a function of lifetime are plotted with sev- 
eral examples of short- and long-lived oscillons, all with initial 
amplitudes Ao = 2. From left to right, the initial radii for the 
oscillons are 2.35, 2.41, 2.53, 2.65, and 2.86. It is quite clear 
that the oscillons decay as they cross the critical values com- 
puted analytically. 



FIG. 10: Oscillon lifetimes vs. decay energy En = Er- Solid 
curve is theoretical, dashed line is numerical. The theoretical 
curve (with an error of ~ 6% in the horizontal positioning of 
the peak) correctly predicts the shape of the distribution and 
that there exists a maximum lifetime in this system on the 
order of ~ 10*. 



Carrying out the above calculations for several trajec- 
tories along the line of existence and plotting the life- 
time vs. Ar, Rr, and Er, results in Figs. [9]and[l0] re- 
spectively. Fig. [9] shows the analytical decay amplitude 
and radius as a function of lifetime (dashed lines) plot- 
ted against several long-lived and short-lived oscillons, 
with excellent agreement: oscillons decay as they cross 
the coordinates specified by the line of existence. Fig. [lO] 
shows the computed lifetime as a function of decay en- 
ergy (E]j = Er), reproducing the maximum lifetime on 
the order of 10'' which is characteristic of oscillons in this 
system. The figure compares the analytical computation 
of lifetime (continuous line) with the numerical results 
(dashed line). The small disparity in the center of the 
peak (of order ~ 6%) is probably due to the Gaussian 
ansatz we use to describe oscillon configurations. 

As was done in Eqs. [34] and [37] our method allows 
us to investigate an oscillon evolving along a particu- 
lar trajectory in detail. As an example, we consider 
the longest-lived oscillon in this model, obtained with 
the initial parameters {Aq — 2; Ro — 2.86). Using 
the information from the curve in Fig. |10[ we can find 
the trajectory whose lifetime corresponds to this oscillon 
(Tiifo ~ 7100), marked as the thicker line in Fig. [S] (the 
coordinates of the reference (or decay) point are Ar ~ .92 
and Rr — 3.25). We then use Eqs. 32 and 33 to calculate 



px and 7x for any parameter of interest. Then, combin- 
ing these values with Eqs. |25] and |D4| we can compute 
the amplitude (Eq. 37 1, radius, frequency, energy (Eq. 
34 1 , and radiation rate as functions of time and compare 
the results with the numerical values. The results are 
plotted in Figs. [TT] [12] and [13] showing excellent agree- 
ment. We can also plot the theoretical prediction for the 




FIG. 11: Comparison of theoretical (continuous line) vs. nu- 
merical (dashed line) radius (top), frequency (middle) and 
amplitude (bottom) for an oscillon with initial conditions 
{Aq — 2, Rq = 2.86), showing very good agreement [theo- 
retical results are computed with {Ar,Rr) — (.92,3.25)]. 



trajectory of this oscillon in the {A, R) plane (shown by 
the thicker line in Fig. Js]), and compare it to the numer- 
ical value, as shown in Fig. [l4] 

In the next section, we will complete the character- 
ization of this system by deriving an expression for the 
frequency of the superimposed oscillation observed in, for 
example. Fig. [9] which seems to be connected with the 
oscillon decay process: the larger the amplitude of the 
superimposed oscillation, the shorter the lifetime. 
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FIG. 12: Comparison of theoretical (continuous line) vs. nu- 
merical (dashed line) values for the energy of the longest-lived 
oscillon, obtained with initial conditions {Aq = 2, J?o = 2.86) 
showing excellent agreement [{Ar,Rr) — (.92,3.25)]. 
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FIG. 13: Comparison of theoretical (continuous line) vs. nu- 
merical (wavy line) radiation rate for an oscillon with initial 
conditions {Aq — 2, Rq — 2.86) showing excellent agreement. 
The inset, which is plotted on a log scale, makes it clear that 
the theory correctly reproduces the rapid initial drop in radi- 
ation rate over many orders of magnitude; the linear scale on 
the larger graph shows that the theory correctly reproduces 
the extremely small (but finite) radiation rate towards the 
end of the oscillon's life [{Ar, Rr) ^ (.92, 3.25)]. 



IV. ANALYSIS OF OSCILLON STABILITY AS A 
FUNCTION OF TIME 

Recall from the first chapter that the three conditions 
representing an oscillon before decay, at the point of de- 
cay, and after decay, respectively, are: 
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FIG. 14: Comparison of theoretical trajectory (solid curve) 
vs. numerical trajectory (dashed curve) for the longest-lived 
oscillon [Aq = 2, Ro = 2.86) showing very good agreement 
during the more stable phase of the oscillon's life {A < 1.5). 
The great increase in density of data points in the dashed 
line in the range .9 < ^4 < 1.5 is due to the prolonged pe- 
riod of time spent in this region by the oscillon (i.e., the 
"plateau" phase). The end point where the dashed line again 
becomes dashed {A ~ .9, 7? ~ 3.2) signals the numerical 
decay point of the oscillon. The thick segment of the solid 
line highlights the portion of the theoretical trajectory dur- 
ing the low-radiation plateau phase; the end of the thick 
segment {Ar,Rr) — (.92,3.24) marks the theoretical decay 
point, showing very good agreement. It is interesting to note 
that, even after the oscillon decay at A ~ .9, the remaining 
field configuration continues to tend to the attractor point at 
{Aoo,Roo) — (.456,4.79), as does the theoretical curve. 



oscillon stability in more depth. In doing so, we will 



obtain a more general formulation of Eqs. 38 which will 
provide a precise measure of the oscillon's stability when 
^gap 7^ riin- We will then see that an expression for the 
frequency of the superimposed oscillation seen in Fig. [9] 
which is clearly related to stability, will naturally emerge. 

In deriving the equations governing the radiation rate 
and lifetime in the previous sections, we made the sim- 
plifying assumption that the oscillons under study are 
long-lived. Mathematically, this assumption is employed 
in approximating the series expansion in Eq. |B7| by its 
first term, yielding Eq. B8 The more stable the oscillon, 
the smaller a given term in the series expansion will be 
relative to the term before it. 

Let iy„ denote the magnitude of the nth term in the 
series of Eq. |B7| In this section, instead of assuming 
that Vn ^ i^n+i, we will compute the fractional difference 
between two adjacent terms and take the result to be a 
natural measure of the stability of the oscillon. 

Define the (dimensionless) stability function E from 



(39) 



(38) two adjacent terms in the series of Eq. |B7| as 



In this section, we seek to investigate the concept of When the oscillon is highly stable, t'n ^ t^n+i, and E 
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1; conversely, as Vn+i ^ i^n (causing the series in Eq. B7 
to fail to converge) then E — > 0. In Appendix E (see Eq 
E7) it is shown that, for any value of n, 



r 



nl 



gap 



(40) 



Eq. 40 is refered to as the stability function; as would 
be expected, it involves the ratio between i^gap and T^i- 
Using Eq. |29]for r„i, we can plot Eq. |40]for the longest- 
lived oscillon (using Eqs. 32 and 33). This is shown in 
the top graph of Fig. 15 The extreme closeness of E 



to unity for most of the oscillon's life, when compared to 
Eq. [39] verifies that we are quite justified in assuming 

The minimum stability allowed at a given time, Emin, 
is attained when T^i is at its maximum, namely, when 

Till — Tlin, 



EmiTi — 1 



r 



lin 



(41) 



Now note that the conditions in Eq. 38 can be written 
in terms of Smin as 




FIG. 15: The top graph shows the theoretical calculation of 
E vs. time and the bottom shows Emin, both for the oscillon 
with {Ar,Rr) — (.92,3.25). The stability measmed by E is 
clearly related to the radiation rate of the oscillon: this kind of 
stability increases in time, since the radiation rate decreases. 
On the other hand, the stability measured by Emin is related 
to the resistance of the oscillon to decay: this kind of stability 
decreases in time as the oscillon moves closer to the line of 
existence. 



Sniin > 0; 
^min — 0, 
Smin < 0, 



(42) 



respectively. In the bottom graph of Fig. [15] we plot 
Emin for the longest hved oscillon (again, using Eqs. 32 
and 33 and with Fiin given by Eq. 



It is clear that one can interpret the stability E as a 
measure of the radiation rate of the oscillon: the radia- 
tion rate decreases in time, so the stability increases. On 
the other hand, Emin measures the resistance of the os- 
cillon against spontaneous decay: as the oscillon evolves 
and approaches the line of existence, this kind of stability 
decreases. 



We will now write the conditions in Eqs. 38 and 42 in 
a third and final form. First, observe that we can write 



as 



where 



"^mod 



F2 



(43) 



(44) 



In terms of Wmodi the conditions in Eq. |42| becomes 

WmodGSR, ^0; (45) 

Wmod — 0; 

respectively. In other words, if we consider the quantity 



decay, at the decay point, and after the decay are 



B{t) = o 
B{t) - 1 
B{t) 



gap 



(46) 



respectively, where the last condition follows since, after 
the decay is initiated, Wgap tends to zero, making c^mod 
tend to zbiFiin. 

Therefore, we can conclude that w^od is a special fre- 
quency associated with the decay of the oscillon whose 
value decreases in time and, at a certain point, be- 
comes imaginary, signaling the oscillon's final demise 
with timescale on the order of the linear decay width 
(Eq. 46 1. In fact, as mentioned previously, such a phe- 
nomenon is commonly observed in, for example. Fig. [9] 
Specifically, there exists a modulation oscillation whose 
frequency tends to decrease as time progresses, until, at 
a certain point, the oscillon decays with width ~ F^n. 

In Fig. |l6]we plot the period 



2tt 



27r 



decay 



^ lir 



(47) 



B{t) 



along with the numerically measured period of the super- 
imposed oscillation. As shown, Tdccay quite accurately 
reproduces this frequency. 

In conclusion, we now have three separate (yet equiv- 
alent) formulations of the condition for oscillon decay. 
The first says that the nonlinear and linear peaks must 
significantly overlap (Eqs. 38 1. The second says that the 



then the real part of B{t) before the measure of oscillon stability must fall to zero (Eqs. 42 1 
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FIG. 16: The solid curve is the theoreticcil CcLlculcttion of ^Tdccay 
for {Ar,Rr) — (.92,3.25). The dashed curve is the numeri- 
cally measured period of the superimposed oscillation, show- 
ing very good agreement. 



FIG. 17: The four oscillon time scales plotted over time. From 

top to bottom we have plotted T^olax, Tdocay, Tlinoar, and Tosc- 

It can clearly be seen here that an oscillon is an object gov- 
erned by multiple timescales spanning many orders of magni- 
tude. 



The last states that the modulation frequency Wn 
become imaginary (Eqs. 451. 



□d must 



V. THE FOUR OSCILLON TIMESCALES 

We will now review the four timescales associated with 
oscillons encountered in our theory. They are: 



'^elax - 
'^ccay 



1 



E -/e[E-E^]3 
2n 



2 

lin 



-'linear ^mass-'*' 

-'CSC 



(48) 

(49) 

(50) 
(51) 



The first, T^oiax, is the relaxation time of the oscillon 
and is typically the longest timescale present. This is the 
timescale over which the oscillon experiences significant 
change. It's net value expresses the inverse rate of energy 
radiation, being thus largest where the oscillon radiates 
the least, as can be seen from the flatness of curves such 
as those in Fig. [TT]and[T2j This is linked to the lifetime 

by, 



life 



7^.9 [^D - E^]i 



(52) 



E—Ec 



[■elax 



E—Ejz 



where we've used Eq. [28] and the fact that the dynamical 
exponents are typically of order unity (see Fig. |6| . 

The second timescale, Tdocay, is the period of the su- 
perimposed oscillations seen in the oscillon as a result 
of its motion towards the line of existence. The third is 



the decay time of an object in the linear theory [S]. The 
fourth, and shortest timescale is the oscillation period of 
the oscillon. 

The theoretical values of the four timescales are plotted 
in Fig. [17] vs. time for the longest-lived oscillon. Note 
how, together, they span many orders of magnitude. It 
is the presence of these four widely different timescales 
in one single system that makes oscillons such intriguing 
objects to study. 



VI. CONCLUSIONS AND OUTLOOK 

In this work, we expanded upon and improved our 
method to compute analytically the decay rate and life- 
time of oscillon- like configurations [IF. Our approach 
relies on the comparison between the radiation spectrum 
of the nonlinear, oscillon-bearing, model and its linear 
limit. The radiation ultimately responsible for the oscil- 
lon decay is related to the overlap between the two spec- 
tra: the larger the overlap, the faster the decay. We have 
shown that in both d = 2 and d = 3 there is an attrac- 
tor point in field configuration space and that oscillons 
may form as configurations evolve toward that point. In 
d = 3, we may think of oscillons as unstable but very 
long-lived perturbations about the attractor point. In 
d = 2, oscillons may be considered as perturbatively sta- 
ble perturbations about the attractor point. All oscillon 
configurations are shown to migrate toward the attractor 
point as they evolve in time. However, in d = 3 this point 
lies below the oscillon's line of existence, and no oscillons 
in the models analyzed here, double- well potentials, are 
absolutely stable. It remains to be seen if it's possible to 
find models in c? = 3 where oscillons are perturbatively 
stable. In c? = 2, the situation is different: the attrac- 
tor point lies within the oscillon region and our theory 
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predict that they should be at least perturbatively sta- 
ble or exceedingly long-lived. We can thus interpret the 
oscillons in symmetric double-well potentials as special 
types of perturbations about the attractor point: those 
in d = 3 are unstable but can be very long-lived, while 
those in d = 2 are perturbatively stable. Comparison 
with numerical results show that our method provides 
a very good quantitative description of these configura- 
tions. We remark that there are very few examples in 
the literature where analytical results provide a good de- 
scription of time-dependent, nonperturbative phenomena 
in relativistic field theories in d > 1. We have also ob- 
tained a precise criterion to establish the stability of these 
configurations, as encapsulated in the oscillon stability 
function of Eq. |4l] 

It should be born in mind that the assumptions used in 
this paper, although quite general, limit the applicability 
of the theory. First, it is known that the oscillon solu- 
tion is not an exact Gaussian. Therefore, given that all 
quantitative predictions made by the theory are based on 
the Gaussian, all results will possess a slight error. The 
exact nature of this error is difficult to quantify; how- 
ever, it should be small given that the oscillon (in the 
models studied here) is known to be well- approximated 
by a Gaussian. Another limitation of the theory lies 
in Eq. D4 In the models studied in this paper, the 
core amplitudes of oscillons with shorter lifetimes oscil- 
late with large amplitude about their mean value. (See, 
e.g., Fig. [9j) Since the energy E{t) does not oscillate, it 
follows that it is not possible to state, for example, that 
A{t) — Aoo — -fA{E — EooY^ , when the oscillation ampli- 
tude of the superimposed oscillon (or of the radius, R{t)) 
is large. Thus, our theory cannot make accurate predic- 
tions of the time-dependence of individual parameters for 
oscillons which are relatively short-lived. 

Given the generality of our method, we expect it to be 
extendable to many other models. For example, a simple 
next step would be to apply it to oscillons in asymmetric 
double well potentials, where their longevity is expected 
to increase [6] . We currently searching for models where 
oscillons in d = 3 may have an attractor point above 
the line of existence. The relative location of the at- 
tractor point with respect to the line of existence should 
serve as a general criterion to determine the longevity 
of time-dependent scalar field configurations in a vari- 
ety of models, including those with more than one scalar 
field. Quite possibly, there may be models with two cou- 
pled fields that produce stable oscillons even in d = 3, 
at least in the sense that Smin ^ 1- Another possible 
extension of the present results is the inclusion of gauge 
fields. As shown in Refs. [TT] and [12], in both U(l) and 
SU(2)XU(1) models oscillons have not been seen to de- 
cay. The question of their absolute stability remains open 
and is of obvious interest. Finally, oscillons may play a 
key role in the dynamics of the early universe: they may 
be formed during preheating after inflation and delay 
thermalization; they may leave behind a gravitational- 
radiation signature; and they may contribute to the dark 



matter component of the cosmic energy density. We are 
currently investigating how to extend the current meth- 
ods to an expanding cosmological background. 
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APPENDIX A: LINEAR RADIATION 
DISTRIBUTION IN d = 2 

We begin by expanding the Gaussian in eigenfunctions 
(zeroth-order bessel function of the first kind) of the two- 
dimensional, spherically symmetric Klein-Gordon equa- 
tion: 

A{t)e-'^ = f b{k)Jo{kp)dk. (Al) 
Jo 

Using the orthogonality relation 



pJo{kp)Jo{k'p)dp = jS{k - k') (A2) 



we can invert Eq. |Al|to yield 



6(fc) = kA{t) / Jo{kp)e-'^pdp. (A3) 

"'0 

In |il7j it is shown that integration gives 



(A4) 



Note that this function has the same form as the corre- 
sponding one in d = 3. Hence the results in Eqs. [6] and 
[7] need not be recomputed. To obtain Fiin in d = 2 one 
could, in principle, solve the equation of motion as was 
done in d = 3; however, for our purposes it will be suf- 
ficient to determine this parameter numerically (by nu- 
merically integrating the linear Klein-Gordon equation in 
d — 2 with a Gaussian initial condition) , yielding 



hn — rj- 
^ ■'line 



/?2 



.6 

W 



(A5) 



APPENDIX B: DERIVATION OF NONLINEAR 
FREQUENCY PEAK 

In this paper we consider spherically symmetric oscil- 
lons which can be accurately modeled by an oscillating 
field configuration whose spatial profile and amplitude of 
oscillation vary little over the course of an oscillation, i.e., 

(/.(r, t) ~ A,{t)P(r- R) = A,,,{t)A{t)Pir; R) (Bl) 
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where (j) is the field, r is radial position, t is time, P{r) 
is the spatial profile of the oscillon normalized so that 
P{r) = 1 at the origin, i? is a time-dependent measure of 
the spatial extent of the oscillon (i.e., the "radius") which 
is assumed to vary little over the period of a single oscil- 
lation, Ac{t) is the time-dependent oscillon core [(/)(0,i)], 
A(t) is its time-dependent envelope of oscillation, and 
^osc(^) is an oscillating function which is normalized to 
an upper turning point of unity. 

If the oscillon oscillates approximately harmonically we 
can write 



^osc(^) 



X _^ COs(cJnl^) 

1+X 1+X 



(B2) 



where Wni is the time-dependent frequency of the oscillon 
and X is a dimensionless constant which accounts for a 
possible non-zero center of oscillation. 

Now assume that we are interested in calculating the 
radiation rate of the oscillon at t = 0. Since the oscillon 
radiation rate depends only on the instantaneous proper- 
ties of the oscillon, writing UJn\{t) = Wni(O) in Eq. B2 will 
still yield the correct radiation rate at t = 0. Combining 
Eqs. |Bl| and [B2] and taking the unitary cosine transform 
of Ar we have 



cos{ujt)Ac{t)dt 

' A(t) 
cos(wt) cos(u;ni(0)) ^ dt. 



(B3) 



1 



X 



where in the second step we've multiplied A(t) by Aosc(i) 
in Eq. |B2| and dropped all terms except the one propor- 
tional to cos{LL!jiit) A{t) which will create a finite width 
peak centered at the oscillon frequency (all others are ir- 
relevant and do not contribute to the radiation rate). We 
can rewrite Eq. |B3|as 



1 



[cos(w-i_t) -I- cos( Lu_t)] ^^^^ dt, 
1 + X 



(B4) 



where w+ = + lu and lo^ = Wni — oj. Since each of the 
two cosine terms above will contribute an identical peak 
(one centered at -|-ajni and the other at — Wni) and both 
will contribute identically to the radiation rate, we can 
simply drop the one centered on —lo-^\ and multiple by 
two, yielding 



cos([lj„i 



j]t)Y^dt. 



X 



(B5) 



Now, since an oscillon's amplitude A{t) tends to a finite 
value, denoted Aoo, as i ^ oo, we define = A{t)—Aoo 
and write 



where we've dropped the transform of the term since 
it will produce a delta function centered at Wni which will 
not in anyway affect J-{uj) in the region of interest (the 
radiation zone). 

To proceed, we use the fact (found be performing suc- 
cessive integration by parts) that, for a function f{x) 
which tends to zero as x — > oo. 



cos{sx)f{x)dx = ^(-l)"^^-li, (B7) 



where /*^™' is the mth derivative of the function /. We 
apply Eq. |B7| to Eq. |B6| and note that, for a long-lived 
oscillon, the flatness of A{t) implies that, for s in the tail, 
we need only keep the first term, yielding 



2 A 



1 + X i^nl - t^)' 



(B8) 



where we've switched the evaluations at t = to a general 
time t, since the above calculation can be applied to any 
physical time. Eq. |B8| essentially states that when a 
peak associated with a decaying function has a very small 
width then the tail goes like (wni — ^)^'^ 



It will be useful to note that ^{10) in Eq. B6 will be a 
peak whose width, denoted Ta, scales as 



A 



A 



X 



d 

di 



Aa 
1 + X 



A 



A- A. 



(B9) 



where we've used that x is a constant, Aa = A — A^ 
and Ana is a constant. 



APPENDIX C: OVERLAP FUNCTION 



Consider the overlap function fl{uj) in Eq. 17 



n{uj) = aT{uj)b{uj) = T{uj)h{uj). 



(CI) 



Since ^{uj) has dim ension of amplitude per unit fre- 
quency (see Eq. 16 1, as does T{uo), it follows that h{uj) 



is dimensionless. Therefore, one can interpret h{uj) as a 
dimensionless coupling factor which modulates the distri- 
bution J-{uj) . If the oscillon couples weakly to a certain 
mode, biuj) will be small at that value of w, etc. 

Now, biuj), which is proportional to the linear radiation 
distribution, will be a peak whose frequency of maximum 
coupling tJniax is found by 



db_ 
du 



0. 



(C2) 



cos([a;ni — i^]t) 



X 



dt (B6) 



cos([a;ni 



uj]t)^^dt, 



Consider now that the largest-amplitude radiation 
wave which can be created by a driving force of fre- 
quency ijj > Wmass will be produced by driving at the 
frequency ti-imax- However, the largest-amplitude radia- 
tion wave which can be created by a given driving force 
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will have an amplitude on the order of the driving force 
itself (never larger). Therefore, at frequencies near Wmax, 
the amplitude of the radiation wave f2 will be roughly 
equal to the amplitude of the "driving force" T . In other 
words, ri(a;max) -^^l^max), implying that 6(wmax) 1 
from Eq. [Cl] This leads to 



1 



a : 



&(^max) 

which, when combined with Eqs. |17| and [B8| yields 



2 A 



1 



1 + X (^nl fe(t^max) ' 



(C3) 



(C4) 



For example, in the system given by Eq. [T| where hiiS) 
is given by Eq. [o] Wmax = ^/2+~2jB? . This leads to 



e^/'^R/\/2 and, when combined with Eq. 



17 



R A (a;2_2)5 

^ — 

v/tt 1 + X {^n\ - ujY 



■('^ -2) 



gives 



(C5) 



This function is plotted in the inset of Fig. [5] for typical 
values of the parameters in the system given by Eq. [ijin 
d — i. Eq. C5 gives the distribution (amplitude per unit 
frequency) of the radiation wave emitted by the oscillon 
in the vicinity of a time t. Its integral with respect to 
frequency gives the total amplitude of the radiation wave 
emitted by the oscillon, denoted A. 

Given an expression for ^(a;), 5 and Wiad can then be 
given by 



f^'(t^rad) = 0; 



(C6) 



/ (CJ - CJrad)^ ri(cj)dw 

J Q.{Lo)duj 



where the prime denotes a derivative with respect to lj. 



APPENDIX D: DYNAMICAL EXPONENTS 

The time dependence of the oscillon and all of the var- 
ious parameters which describe it is due to a single phys- 
ical mechanism, namely the emission of radiation. If the 
radiation rate is high (low), all parameters will change 
rapidly (slowly). In other words, the general nature of 
the time dependence of a parameter will follow that of 
any other parameter. 

In analogy with Eq. 



B9 define 



X 



(Dl) 



which can be interpreted as the decay width for the pa- 
rameter X. Now consider the generic parameters X{t) 
and Y{t) (which could represent amplitude, radius, fre- 
quency, etc.). Intuitively, we can say that the timescales 
Fx and Vy associated with the rates of change of the 



parameters X{t) and Y{t), respectively, are to be pro- 
portional. 



Tx oc F-y 



(D2) 



for any parameters X{t), Y(t) (i.e., if X{t) slows down 
by a factor of two, then so will every other parameter). 
We can then write, quite generally. Fx esc F^; for any 
parameter X, where E denotes energy. Combining this 
with Eq. Dl yields 



X 



E 



X{t) - X, 



Px 



E{t) - E^ 



(D3) 



where px denotes the proportionality constant. Integrat- 
ing both sides of Eq. |D3| yields 

(D4) 



[Xit)-X^]^jx[Eit)-E^]P-, 



where 7x is the constant of integration. 

We will now derive an expression for the general expo- 
nent Px- First define the new variable 



iyx=HX~X^), 



(D5) 



which wc will henceforth employ in this Appendix, in- 
stead oi X, to describe the oscillon. When combined 
with Eq. Ids] this gives 



px 



dvx 

dvE ' 



(D6) 



During the evolution of an oscillon, the change in the 
coordinate i^x is given by 



dux = Vj^x ■ dv, 



(D7) 



where dv is a differential vector which lies tangent to 
the trajectory of the oscillon in i7 space and Vi^x is the 
vector gradient of vx whose direction lies perpendicular 
to the "level" curves associated with z/x- To calculate 
dv we note that, by virtue of the attractor-like nature of 
oscillons, the oscillon will evolve according to a trajectory 
which runs perpendicular to the level curves associated 
with the radiation rate. Mathematically, we can write 



Vi'^ X dp = 0. 



(D8) 



To proceed, we will make use of our initial assumption 
in Eq. Bl that the oscillon (and hence any oscillon pa- 



rameter X or vx ) can be taken to be a function of two 
degrees of freedom. These two degrees of freedom can 
be arbitrarily chosen to be any two independent oscillon 
parameters. Choosing the coordinate pair i> = {vet^a) 
Eq. 



D8 becomes 



OVE 

Combining this with Eq. 
have 



p^duE ^ 

OVA 



(D9) 



D7| and dividing by dvE we 



dvx 
duf 



dvx 

dvE dvE 



dt^E dvx 



duE 



(DIO) 
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Combining Eqs. 24 and D6 we have 



From Eq. B7 we have, 



p^ = 2- 2pA + Pr) 

d 



(Dll) 



- — {2vE ~ 2vA + v,^) = - — v^, 
di>E dvE 



which impHes that, up to a constant, Vj^ = ivE^^VA+v-q- 
Substitution of this expression for into Eq. DIO and 



combining with Eq. D6 yields the desired result: 
(2+fe)lS + (fe-2) 



Px 



dux 

duA 



(2+ Is) 



(D12) 



APPENDIX E: DERIVATION OF STABILITY 
FUNCTION 



From Eq. |D4| it can be shown by differentiating that, 

X^{X-X^Y (El) 

where a = 1 + By further differentiating, it can be 
shown that 



X 



\X-X^ 



= r^-", (E2) 



where X^™) denotes the mth derivative of X with respect 
to time. Now, from Eq. |39l 



(E3) 



I'ri+l /(2"+l) / /(2n-l) /A(2»-l) 



(E4) 



where we've used Eq. B8 and the fact that x is a con- 
stant. Combining Eq. E4 with Eq. E2 (in the case that 
X ^ A) we have 



^^e^^'^^f^^y, (E5) 



\ J \ Fnl 

where /3 is a proportionality constant and we've used that 



Fnl = F^. Combining Eq. E5 with Eq. E3 we have 



E = 1 



r„i 



1 



P \UJ- w„i 



(E6) 



Eq. E6 is a function of uj; to determine the appropriate 
value of Lu we note that, if F^i were to attain its maximum 
value of Fiin and if Fun = Wgap, the oscillon would decay; 



therefore, S = when Fni 



lin 



This implies 



that P{uj — LUni) = ^gap- Substitution of this into Eq. 
|E6| yields the desired result of 



nl 



gap 



(E7) 
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